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Introduction: Pahoehoe lavas are recognized as 
an important landform on Earth, Mars and Io [1,2]. 
Observations of such flows on Earth (e.g., Figure 1) 
indicate that the emplacement process is dominated by 
random effects. Existing models for lobate a' a lava 
flows (e.g., [3-6]) that assume viscous fluid flow on an 
inclined plane are not appropriate for dealing with the 
numerous random factors present in pahoehoe em- 
placement. Thus, interpretation of emplacement condi- 
tions for pahoehoe lava flows on Mars requires funda- 
mentally different models. A new model [7] that im- 
plements a simulation approach has recently been de- 
veloped that allows exploration of a variety of key 
influences on pahoehoe lobe emplacement (e.g., source 
shape, confinement, slope). One important factor that 
has an impact on the final topographic shape and mor- 
phology of a pahoehoe lobe is the volumetric flow rate 
of lava, where cooling of lava on the lobe surface in- 
fluences the likelihood of subsequent breakouts. 

Basic Approach: To simulate pahoehoe lava em- 
placement, Glaze and Baloga [7] considered the 
movement of small parcels of lava with a volume equal 
to the size of a typical toe (70 x 70 x 20 cm 3 ) [8]. The 
model provides a set of probabilistic rules for deter- 
mining the location and direction of movement for 
each parcel. For a constant volume flow rate feeding 
the lobe, a single parcel is added at each time step. 
However, unlike the classical random walk of Browni- 
an motion [9], only one parcel is allowed to move at 
each time step. Thus many parcels may remain 
dormant, but fluid, for multiple time steps. The net 
effect of this new type of random walk is that parcels 
tend to accumulate preferentially within the lobe pro- 
ducing cross-sectional topographic profiles with a me- 
dial ridge (Figure 1), very similar to inflated pahoehoe 
lobes observed in the field [10]. Because lava transfers 
within the lobe dominate the probabilistic activity, this 
approach provides a basis for simulating pahoehoe 
inflation in terrestrial and planetary settings. 

Influence of Volume Flow Rate: The objective of 
this current study is to characterize how volume flow 
rate influences the shape (both planform and cross- 
sectional topography) of a pahoehoe lobe. In the origi- 
nal algorithm [7], all positions occupied by lava at a 
particular time step were considered equiprobable lo- 
cations for the subsequent parcel transfer. The basic 
framework is based solely on stochastic behaviors and 
contains no physical parameters other than flow rate 
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Figure 1. Example of a purely random reference simula- 
tion case with 500 parcels released from a point source. 

(a) plan-view of lobe (enlarged and smoothed for graph- 
ical interpretation), (b) topographic profiles taken 
through the origin along the x- andy-axes. 

and parcel size. The new approach described here in- 
jects the process of surface cooling and skin growth 
into the model by relating measured or computed ex- 
pressions for cooling into the probability of a breakout. 

Concept of Correlation. “Correlation” is a statisti- 
cal term that describes the influence of prior steps. 
When correlation is present, each lava parcel transfer 
retains some memory of earlier time steps and is no 
longer completely independent of prior transfers. Glaze 
and Baloga [7] used correlation to describe sequential 
breakouts at the margin of a pahoehoe lobe. Based on 
the logic that a recently active parcel is more likely to 
be the source of the next breakout at the margin, the 
model arbitrarily increased the probability of a sequen- 
tial breakout (additional parcel added at the same loca- 
tion in the same direction) for the most recent parcel 
transfer [7]. The incorporation of correlation in this 


way produces simulated lobe morphologies that are 
consistent with data collected in the field. 

Cooling Curves and Probability Distributions. De- 
spite the encouraging results of including correltation 
in the random simulation model, the probability distri- 
bution applied to the selection process by [7] was not 
based on any fundamental physical constraints. Here 
we emphasize the role of skin cooling as the dominant 
physical process determining the shape of a pahoehoe 
lobe. The basic assumption of this approach is that 
once a cell initially becomes occupied, the parcel be- 
comes fixed on the surface of the lobe. The mechani- 
cal strength of the crust grows with time and is thus 
proportional in some sense to cooling. Crisp and 
Baloga [11] and many others derive the skin depth, 
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that gives the proportionality. An increase in the resi- 
dence time on the surface decreases the probability of a 
breakout at that location due to the growing thickness 
of the crust. Thus there is a relationship between 
measured cooling and the probability of a breakout. 

The cooling curve shown in Figure 2 is thus direct- 
ly related to the probability of a breakout when proper- 
ly calibrated by the flow rate and normalized. The crust 
is no doubt significantly thinner for t 8 and t 7 than for t 2 
and ti so they are more likely to be the sites for the 
next breakout. The curve with open diamonds in Fig- 
ure 3 shows the distribution of probabilities for deter- 
mining the location of the next time step in the series 
when cooling is considered (t 9 for the series shown in 
Figure 2). The flat curve in Figure 3 (blue diamonds) 
illustrates the equiprobable case for comparison. It can 
easily be seen in the cooling case that the probability 
of the next parcel transfer originating from parcel t 8 is 
greater than from fi. This approach is broadly applied 
to larger simulations to explore the effect of correlation 
on the planform shape and topography. 

Conclusions: Existing models for lava flow em- 
placement based on viscous flow on an incined plane 
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Figure 2. Based on Figure 1 from [12] illustrating field 
observations of surface cooling. Subscripted t’s indicate 
the time step when each parcel (1-8) was effused from 
the source for a particular constant volume flow rate. 
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Figure 3. Probability distribution prescribed for each 
occupied cell location after 8 time steps. 


are inappropriate for application to pahoehoe lavas 
where random processes are dominant. A completely 
new simulation approach developed recently provides 
a framework for exploring the effects of cooling on the 
final morphology of pahoehoe lobes. The basic ap- 
proach assumes that the likelihood of parcel transfers 
is correlated with the surface temperature, with warmer 
parcels (most recently emplaced) having higher proba- 
bilty and older parcels being less likely. 

As suggested by Figure 3, the influence of cooling 
on breakout locations decays with time as more cells 
and potential locations are added to the lobe, i.e., the 
fraction of hotter to cooler parcels decreases with time. 
Unless other factors cause a correlation in the lava 
transfers, the topographic profiles will always ap- 
proach those of the uncorrelated cases described in [7]. 

Results to date indicate that a great deal of infor- 
mation about the emplacement of inflating pahoehoe 
lobes can be extracted by stochastic modeling of the 
random effects. Refinement of this approach will have 
implications for the analysis of inflated pahoehoe lobes 
on the Earth, Mars, Io, and possibly the Moon. 
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